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1. INTkODUCTION 


The Profile Similarity technique for crop classification developed by 
Dr. Gautam Badhwar of the National Aeronautics at«d Space Administration 
(NASA), Lyndon B. Johnson Space Center (JSC), has proven effective for the 
classification of corn (ref. 1). This method Incorporates the effects of 
emergence date distribution Into the classification and bases classification 
on the temporal profile of the crop of Interest (refs. 2 and 3). The 
classification method Is for a specific crop. ResulU of applying the 
technique to the classification of spring wheat In the U.S. northern Great 
Plains are documented In this report. 

The procedure used to apply software programs developed for the classification 
of com to the classification of spring wheat Is given, and numerical results 
are presented. The site data set Is listed In section 2, along with the 
Accuracy Assessment (AA) ground- truth percentages for spring wheat, barley, and 
oats for each site. Section 3 describes the procedure followed for segment 
classification. Results of classification. In tabular form, are presented In 
section 4. Concluding remarks are given In section 5, and a recommended 
procedure for operational use of this technique Is presented In section 6. 


2. DATA SET 


The segment data set for this study consists of 17 spring wheat sites for 
which the AA digitized ground- truth maps are available. This data set was 
used for the Analyst Labeling/Procedure M experiment (ref. 4).^ The set has 
also been processed using Procedure 1 (ref. 5). 

reographical distribution of the sites Is limited: 13 sites are In North 
Dakota. 3 are In Minnesota, and 1 Is In Nebraska. Crop year 1978 acquisitions 
from landsat-2 and Lan(isat-3 were used. The full segment, 22 932 picture 
elements (pixels), was classified. 

Table 2-1 lists the sample segment nunbers for the sites; the location of the 
sites; and the AA ground-truth percentages of spring wheat (Including durum 
wheat), barley, and oats In each of the segments. 


^ 

The origins* data set contained 18 sites; however, becbuse the AA digitized 
grounc*. truth was not available for segment 1835 (Ottertall, Minnesota), this 
site was omitted from classification. 
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TABLE 2-1.- AA GROUNO-TRUTH PERCENTAGES FOR THE 
SAMPLE SEGMENTS USED IN aASSIFICATION 


Sample 

segment 

Location 
(county, state) 

Spring wheat 
(Including 
dunm wheat) 

Barlpy 

Oats 

1380 

Kimball. Nebr. 

7.00 

|KQ|B 

1.81 

1392 

Benson. N. Dak. 

26.24 

5.37 

1.14 

1457 

Ward.'N. Dak. 

42.04 

?.22 

2.67 

1461 

Pierce. N. Dak. 

31.70 

4,67 

3.47 

1467 

Towner. N. Dak. 

39.81 

10.79 

0.31 

1473 

Cass. N. Dak. 

31.73 

16.99 

0.64 

1518 

Roseau. Minn. 

22.19 

2.79 

7.53 

1566 

Red Lake. Minn. 

17.70 

5.19 

5.07 

1602 

Mountrail. N. Dak. 

26.47 

1.08 

1.90 

1612 

McHenry. N. Dak. 

10.99 

0.26 

0.23 

1619 

Grand Forks. N. Dak. 

35.72 

0.41 

0.30 

1636 

Stutsman. N. Dak. 

36.76 

2.24 

3.90 

1650 

Hettinger. N. Oak. 

16.30 

0.91 

4.43 

1653 

Burleigh, N. Dak. 

14.64 

0.40 

3.71 

1656 

Morton, N. Dak. 

3.75 

0.47 

2.85 

1825 

Norman, Minn. 

12.83 

4.88 

8.45 

1920 

Sioux, N. Dak. 

16.89 

0.47 

4.90 
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3. PROCEOUKE 


A detailed description of ttie modeling used In Badhwar Profile Similarity 
classification is given In references 2 and 3. Implementation methods are 
essentially the same as those explained In reference 1 for corn. All data 
processing used to generate the classifications was done on the progranmed 
data processor. Model 11/45 (POP 11/45). The software programs which are 
referred to in the procedure below are described in appendix A of this 
document. 

a. Large Area Crop Inventory ixperlment (LACIE) segment Images for all avail- 
able acquisitions were unloaded to a POP 11/45 disk using the IMUNLD2A 
program. 

b. The data quality of the Imagery for all available acquisitions was noted 
using the production film converter (PFC) film products of these 
acquisitions. 

c. Using the PFC products, at least four candidate training fields of spring 
wheat were defined with a tentative set of acquisitions. Reference to 
ground truth confirmed the field labels. If possible, AA special fields 
were Included. 

d. Candidate training fields were graphed over the tentative acquisition set 
using the IMAPLT software program. If necessary, the fteld definition was 
revised. Occasionally, signature abnormalities which coaid not be seen on 
the Imagery eliminated a field from use, and new fields were defined. The 
acquisition choice was verified from the IMAPLT graphs; If necessary, 
alternate acquisitions were selected. 

e. Two fields were selected as training fields; these fields were those which, 
based on the Imagery and graphs, could be expected to produce acceptable 
classifications. The remaining fields were used as test fields. Since tr.e 
study was conducted as a research and development effort, this field selec- 
tion step sometimes Involved several Iterations using the IMAPLT program. 
The sites which were used exhibited a range of problems normal for LACIE 
segments In the U.S. northern Great Plains; e.g., low percentages of wheat. 
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strip fields* Inadequate and poorly distributed acquisitions* and poor data 
quality. The effects of these problems on classification Mere Informally 
assessed; none precVjded the generation of acceptable classification 
results. HoMever* defining suitable training data* Mhich Is essential for 
good classification* was made more difficult by the problems evident In 
one-half of the sites used for this study. 

f. For each of the training fields* a data file was established which 
contained the acquisition set* training field coordinates (line* pixel)* 
ana initial guesses (derived from the IMAPLT graphs) for the modeling 
constants for the training fl!*ld. 

g. Classification was done In the batch-processing mode. The required 
computer time varied from 20 to 60 minutes. 

h. From the classification flles^ generated based on each training field* 
statistical summary sheets were output. The map formats of the classifi- 
cation file and the digitized ground-truth file were compared. These 
files were evaluated as described In reference 1. Areas of disagreement 
were examined; and. If possible* the reason for the disagreement was 
Identified. 

I. The test field classification was noted. This proved to be useful for 
assessing the classification results. 

J. Classification results were compiled Into the tables presented In 
section 4. Detailed records of the processing done for each segment are 
available. 


9 " — ■ ■ ■■ 

‘’The classification files generated for this report have been released to AA 
for comparison with results using the Procedure 1 classification technique. 
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4. CLASSIFICATION RESULTS 


Classification results are presented In tabular fom. luo classifications, 
based on different training fields. Mere produced for each segoent. For sone 
segments, ttw classification results are very different. These differences 
may be due to the extreme shortage of suitable training fields or to data 
quality problems. Field 1 Is the preferred training choice for the segment. 
Results are given for AA pure pixels only and for all (22 932) pixels.^ 

All percentages are as calculated by the NISMAP program. The percentage of 
the segment not Identified by ground truth plus the classified and rejected 
percentages will be 100 percent of the segment. 

For this study, spring wheat Is defined as spring wheat and durum wheat only 
(AA codes 95. 100. 120. and 125). Other small grains are considered 
misclassifled If they classify as spring wheat. 

For each segment, the following Items are listed In table 4>1. 

a. Sample segment number 

b. Sample segment location 

c. Acquisitions available (Julian date) for the segment [Consecutive-day 
acquisitions are omitted. Landsat-3 acquisitions are denoted by (3).] 

d. Acquisitions used for the classification results presented 

e. Coordinates of each of the fields used to train the classifier 

f. The nunber of pixels In each training field 


the percentages based on pure pixels (defined by AA to be those pixels which 
on a subpixel level contain only oim crop) appear above the dash^ line In 
table 4-1; for all pixels, below the dashed line. 
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g. A confusion matrix of the classification in scene percentages: 

(1) ground-truth spring wheat classified as spring wheat (S ♦ S), 

(2) ground-truth spring wheat not classified as spring wheat (S ♦ N), 

(3) ground-truth nonspring wheat classified as spring wheat (N ♦ S), and 

(4) ground-truth nonspring wheat not classified as spring wheat (N ♦ N) 

h. Proportion estimate comparison of the classified proportion of spring 
wheat and the ground- truth proportion (These proportions do not include 
the area which was not identified by ground truth.) 

i. The percentage of the segment not Identified by ground truth [This 
includes unknown fields (AA code 80), unidentified areas (i^A code 164), 
and areas block-identified as str^p fields of spring wheat (AA codes 170, 
175, 220, and 225).] 

j. Additional comments, including the number of AA pure pixels in the scene 
which was used to compute the percentages 
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TABLE 4-1.- CLASSIFICATION RESULTS 
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TABLE 4-1.- Continued 
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S. CONaUSIONS 


The average percentages of mlsclasslflcatlon In the ground-truth-identified 
area of each segment using field 1 are 17.6 percent for pure pixels only and 
19.8 percent for all pixels. Using field 2 , misclasslflcatlons average 
25.6 percent for pure pixels and 27.4 percent for all pixels. Matrices of 
average mlsclasslflcatlon are shoun below. 


Training 

field 

Pure pixels 

All pixels 

1 

mm 

8.7\ 

/12.4 

10.4 \ 


IQ 

65.1/ 

\ 9.4 

62.0/ 

2 

wM 

8.4\ 

/12.9 



IIH 

56.8/ 

\l7.4 

54.0 / 


A scene accuracy of 75 to 80 percent Is a reasonable expectation for the 
classification of spring wheat using the Badhwar Profile Similarity technique. 

This method of classification can be applied effectively to segments with a 
very low percentage of wheat. Finding suitable training fields In these 
segments can be difficult; however. If a good crop profile Is defined, 
classification results are good. 

The chi-square value gives an estimation of the adequacy of the crop profile 
curve as an approximation of the training field data. It must be used with 
visual examination of the training field data to assure that (1) a curve is 
defined and (2) the field data are compact. Chi-square values are dependent 
upon the standard deviations of the data and must be monitored. This Is 
cumbersome. More objective criteria should be provided for assessing train- 
ing field data approximation. 

The use of test fields Is an efficient and effective aid to assessing 
classification results. 
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For North Dakota* the acquisition coverage and the range alloued for estimated 
planting date In the program (currently tlS days from the estimated planting 
date of the training field) should be extended. In North Dakota* spring iiheat 
has an early and a late planting. If acquisition coverage of the segment Is 
cut off In August* after the early planted Mheat Is harvested* the late 
planted wheat tends to be misidentifled as a summer crop. The option to 
remove the restriction on the planting date range should be provided. 
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6. KECOMCNDATION 


Spring wheat In these segments In the U.S. northern Great Plains was clas- 
sified using titt same program ~ with no changes — used for the classification 
of com In the Com Belt*^ The overall accuracy of the results confirms the 
adaptability of the Badhwar Profile Similarity classification technique to a 
variety of crops. This accuracy also Indicates that the software program used 
for classification, CLASFYT, should be considered operational for spring 
wheat, as well as for corn. 

The following procedure is recommended for operational use of this 
classification technique. 

a. Choice of training field and test fields: 

VIheat Is distinguished from other vegetation by Its growth cycle over a 
time Interval*. In the Landsat bands covering the visible spectral 
regions, channels 1 and 2, this cycle defines a curve similar to that 
shown below In reflectance versus time. 

o 
c 

5 

«*- 

6 

Time 

In the near Infrared spectral regions, channels 3 and 4, the curve Is 
similar to the following: 

c 

& 

Time 

^Thls work Is being documented by the author at the present time. 
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On the PFC Imagery* where channeU 1 and 2 are blue and green and channel 4 
Is red, this relative channel reflectance change depicts a time signature 
change of gray (bare soil; approximately equal spectral response and color 
components) to red (vegetation; since the Infrared value rises as the 
visible value decreases) and back to yellow or gray tones (stubble and bare 
soil). An appropriate time Interval can be determined from ancillary data; 
I.e., crop calendars and regional statistics. 

e Choose four candidate fields which are (1) of a 20-p1xel minimum size; 
(2) free of roads, nonagricultural components, etc.; and (3) free of 
clouds and haze. These fields should exhibit a continuous gray to red 
to brown/gray/green signature sequence In the proper time Interval. 

Vw PIC produet 3 . 

e Determine field borders; border and edge pixels must be avoided In 
field definition. Define field coordinates at least two pixels Inside 
field borders on all available acquisitions. Ueo PFC produete 1 and 2, 

e Designate one field as the training field on which the classification 
will be based. The remaining fields will be used as test fields to aid 
In classification evaluation. 

b. Choice of acquisition set ; 

Acquisitions for classification are chosen to characterize the wheat growth 
cycle; the set selected should be well distributed over this cycle. Cloudy 
or hazy acquisitions, as well as those which are preemergent for wheat or 
which exhibit appreciable amounts of wheat harvest, should be avoided. 
Classification can be done on a set of four or five acquisitions. A five- 
acquisition classification Increases crop separability, but this many 
suitable acquisitions may not be available. 

If problems with the acquisition choice or with the field definition exist, 
the analyst may choose to plot the field and acquisition set using IMAPLT 
before entering the classification. 
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c* Clastifiottlofit 


Classification will ba dona In tha batch-processing mode on the Interactive 
Nultlspectral Inage Analysis System* Nodal 100 (Image 100)* POP 11/45 Image 
processor. The analyst iflll Input the field coordinates and acquisition 
set via cards. The required computer time Is 20 to 60 minutes. 

d. Evaluation of results: 

Products of classification will be: 

e A summary sheet of Input values and calculated parameters based on the 
training field data with a numerical sunmary of'pllcel classification 
results (fig. 6-1) 

a A full-scene classification map generated as a film product for analyst 
use (fig. 6-2) 

The film product classification map should have evident field patterns. 
Fields should be well filled out with a minimum of blank spaces (pixels 
rejected as spring wheat) In the field Interiors. Blank areas also should 
be clear; I.e., be reasonably free of scattered pixels classified as 
spring wheat. Scattered pixels or blank areas may Indicate an overclassl- 
flcatlon or an underclassifl cation. The classification map should be 
overlaid on the PFC film products to check agreement of the classification 
with the analyst's Identification. 

On the summary sheet, calculated parameters should be checked; the esti- 
mated planting date of the training field as generated on each channel 
should be the same within the estimated planting date error, and the 
chi-square fit should be less than 10 In each channel.^ A scene 


^re objective criteria for evaluating the statistical output need to be 
determined. Currently, visual examination of the channel graphs Is used to 
assess the compactness of the training field data and the adequacy of the 
curve approximation of the data. With this, a chi-square fit of less than 10 
Is meaningful. The chi-square fit value Is dependent upon the data variation; 
an evaluation value which combines these parameters should be defined for 
operational use. 
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proportion ostlinato can be computed from the pixels classified as spring 
idieat» and this estimate can be compared with available statistics. As an 
addition to the summary sheet of figure 6-l» the numerical proportion of 
pixels correctly classified as spring wheat should be supplied for each of 
the analyst-defined test fields; the test field accuracy should be 
70 percent or better. 

If the classification Is unsatisfactory* rework will consist of selecting 
an alternate training field or acquisition set. This choice should be 
graphed, using the IMAPLT program, before use as the basis for segment 
reclassification. 
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ACCURACY ASSESSMENT QUALITY ASSURANCE CLASSIFICATION 

PRfCESSING DATE - 12-NAR-80 AT 02:58:04 

SEGMENT NUMBER - 1653 CRIP OF INTEREST - SWHE 

IMAGE FILLS USED IN CLASSIFICATION - DB2:tlll.3ll65378136.IM2 

082: 111,3 165378154.IM2 
082: 111,3 165378191. IM2 
082: 111,3 165378191.JM2 
0B2:ail.3]165378208.IM2 

TRAINING FIELD - LINE N|. SAMP. N|. 

73.0 112.0 

70.0 129.0 

• 73.0 130.0 

75.0 112.0 

MEANS AND STD. DEV. F|R TRAINING FIELD BASED IN 40 PIXELS • 
CHANNEL ACQUISITION OATES 


NUMBER 


78136 

78154 

78191 

78191 

78208 



1 

MEAN 

STD. DEV. 

30.11 

2.17 

23.03 

1.51 

21.17 

0.93 

21.17 

17.76 

23.21 

0.99 



2 

MEAN 

STD. DEV. 

33.46 

2.32 

17.33 

1.73 

18.25 

1.97 

18.25 

37.74 

24.71 

1.63 



3 

MEAN 

STD. DEV. 

37.36 

2.01 

40.83 

2.18 

45.82 

1.74 

45.82 

35.88 

40.61 

2.26 



4 

MEAN 

STD. DEV. 

31.63 

1.63 

36.80 

1.09 

44.35 

2.98 

44.35 

61.04 

36.88 

1.91 



CONSTANTS FIR N|Da • 
CHANNEL 

NUMBER A 

ALPHA 

BETA 

T| 


CHISQ 

I INITIAL 
FINAL 

3.80 

3.77 ♦. 5.75 

.5.37 

-5.51 

♦- 1.65 

-0.82 

-0.84 

- 0.28 

1.20 
1.19 ♦- 

2.19 

0.11 

2 INITIAL 
FINAL 

5.00 

4.69 8.94 

-15.19 

-15.16 

♦- 2.64 

-2.49 

-2.49 ♦- 0.45 

1.20 
1.15 +- 

1.21 

0.97 

3 INITIAL 
FINAL 

3.22 

3.31 ♦- 5.40 

3.91 

3.47 

♦- 1.52 

0.62 
0.54 ♦ 

- 0.26 

1.20 
1.15 ♦- 

3.05 

1.24 

4 INITIAL 
FINAL 

3.00 

3.09 ♦- 0.00 

5,23 

4.17 

♦- 0.00 

0.81 

0.63 0.00 

1.20 

1.16 4-.0.00 

2.22 

CHISQ TWIESHILO - CHANNEL 

1 

2 

3 

4 





TWESHUD 7.04 7.04 7.04 8.17 


CLASSIFICATIIN RESULTS - 

PIXELS CLASSIFIED SW« - 9689 CUT BY CW - 4456 
PIXELS SCREENED - 0 CUT BY CK3 - 3521 
PIXELS CLASSIFIED N|N.SUHE - 13243 CUT BY CHI - 5266 


Figure 6>l«> Example of a statistics summary sheet. 
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Figure 6-2.- Filin product classification map of the full scene. 
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APPENDIX A 

SOFTWARE PROGRAMS USED FOR CLASSIFICATION OF SPRING WHEAT 
USING BADHWAR PROFILE SIMILARITY TECHNIQUE 


APPENDIX A 

SOFTWARE PROGRANS USED FOR CLASSIFICATION OF SPRING WHEAT 
USING BADHWAR PROFILE SINILARITY TECHNIQUE 


hA IHUNLD2A 

IMUNL02A takttt an Ifnaft unload tapt ganerated on tha Earth Rasourcat Intar- 
actlva Procaaslng Sjystan (ERIPS)» adits It usInQ SCREEN (raf. 6)» adjusts tha 
Landsat-3 acquisitions Into a data ranqa comparabla to tha data ranqa of 
Landsat-2 acquisitions using the Wahsianan aultlpllcatlva factors (raf. 7)i and 
loads tha Images on a POP 11/45 d1sk< 

a. Input: ERIPS Image urload tape 

b. Output: scraanad Images with adjusted Landsat-3 acquisitions on a 
POP 11/45 disk 

A.2 IMAPLT 

IMAPLT (ref. 8)^ plots the Individual pixels of a field* giving reflectance 
values versus time (I.e.* the acquisition dates specified) for each channel. 
IMAPLT then plots the field mean values* each channel* with a one standard 
deviation envelope; a curve Is fitted through the mean values. Eight graphs 
(two for each Landsat channel) are produced for a field over a set of acquisi- 
tions. Graphs are displayed on the Image 100 Tektronix screen* and hardcopies 
are made automatically. The segment minber* the acquisitions used* the coor- 
dinates of the field* the channel number* the numb^Nr of pixels In the field* 
and the mean and standard deviation on each acquisition are listed on the first 
plot. The constant values computed from the data for the model (with the 
estimated error)* the estimated planting date of tha field (with error)* the 
values of the fitted curves at the specified acquisitions (which can be 
compared with the computed mean values of the data)* and the chi-square value 
for the fit of the approximating curve to the field data are presented on the 
second plot. 

^Available reference Is to TRJPLT* an ea*^1y version of IMAPLT. 
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a. Input: field coordinates In order; acquisition set of four or five 
acquisitions 

b. Output: . ]ht graphs as described above 
A.3 CLASFYT 

CLASFYT (ref. 9) computes the constants for the curves from the training field 
data, compares the values for each pixel In the segment with these curves (the 
crop profile In each channel),^ and rejects those pixels which are not within a 
specified chi-square measure of the profile. The technique for rejection is to 
compare the pixel channel values with the profiles In channel 2, channel 3, 
channel 4, then channel 1 in succession and reject the pixel If the comparison 
In any single channel is unsatisfactory. Variability of the time of planting 
and/or emergence Is allowed for in the comparison of individual pixels with the 
crop profile (refs. 2 and 3). Accepted pixels are labeled as spring wheat; 
rejected pixels, nonspring wheat. 

a. Input: four or five image files, coordinates of one crop-of- interest field 
to establish crop profiles, and initial values for the function constants 
as computed In IMAPLT (to aid convergence of the approximating curve) 

b. Output: classification file on disk which has a designation of spring 
wheat or nonspring wheat for each pixel in the segment; line-printer sheet 
summarizing the following; 

e Acquisitions used 

e Training field coordinates and the number of pixels in the field 

t Mean and standard deviation for each channel and each acquisition 
(field averages) 

t The input and the final constants (with error) for the model 


^As each Image was unloaded from an ERIPS image unload tape onto a disk for 
processing on the PDP 11/45, It was edited using the ERIM program SCREEN. 
Pixels in the training field failing to pass this edit step were excluded 
from processing, hence from affecting the crop profiles. However, screened 
pixels were restored before classification of the segment so all 22 932 
pixels are designated "wheat" or "nonwheat" (ref. 10). 
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• Final chi-square values for each channel (training field data) 

• Estimated planting date of the training field (with error) as derived 
for each channel 

e Chi-square thresholds In each channel applied as cutoff values In 
classification 

e The nunber of pixels cut for exceeding the chi-square threshold, hence 
removed from consideration as spring wheat. In each channel 

e The final numerical results: the nunber of pixels classified as spring 
wheat, the number of pixels screened (always zero In this study), and 
the number of pixels rejected as spring wheat 

A.4 A2S6MAP 

A2SGMAP provides a full-scene classification map (22 932 pixels) of the results 
obtained using CLASFYT. The scale Is the same as that used for the AA 
digitized ground-truth maps. Pixels classified as spring wheat are designated 
"C" (crop of interest), and those rejected as spring wheat are left as blank 
spaces on the map. 

a. Input: classification file from CLASFYT 

b. Output: line-printer map of the full-scene classification 

A.5 TAPEOUT 

TAPEOUT (ref. 11) reads the data files produced by CLASFYT and creates 
Universal-formatted tapes. Black-and-white film product classification maps 
are produced on the PFC from these tapes. The scale used Is the same as that 
of the PFC color imagery. 

a. Input: classification file from CLASFYT 

b. Output: black-and-white classification map of the full scene on film 
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A.6 MI SWAP 


MI5MAP (ref. 9) compares the classification file produced by CLASFYT with the 

AA digitized ground-truth inventory map for the segment. A line-printer map 

with the following codes is generated: 

a. Ground-truth spring wheat classified as spring wheat appears as S. 

b. Ground-truth nonspring wheat rejected as spring wheat is left blank. 

c. Ground-truth nonspring wheat classified as spring wheat appears as 

d. Ground-truth spring wheat rejected as spring wheat appears as -• 

e. Pixels for which ground truth is not available but which are classified as 
spring wheat appear as $• 

f. Pixels for which ground truth is not available but which are rejected as 
spring wheat appear as %. 

A numerical scene summary Is given in confusion matrix form. 

MISMAP line-printer maps can be generated for all pixels or for AA pure pixels 

only. 

a. Input: classification file from CLASFYT and ground-truth inventory map 
file 

b. Output: full -scene line-printer map comparing the classification file 
produced by CLASFYT with the AA digitized ground-truth Inventory map and a 
confusion matrix nunerical summary of results 
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